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Introduction: 

A  leading  cause  of  cancer  deaths  among  women  worldwide  is  due  to  breast 
cancer  and  new  therapies  and  optimal  treatments  are  continuously  being  conceived  and 
explored  to  better  control  or  even  cure  this  disease.  (1).  Diagnosis  and  therapeutic 
management  of  the  breast  tumor  remain  significant  medical  challenges,  hence  early 
detection,  diagnosis,  and  timely  treatments  are  essential  to  successful  health  care  (2). 
Magnetic  resonance  imaging  and  spectroscopy  (MRI/MRS)  have  gained  significant 
importance  during  the  last  fifteen  years  for  the  diagnosis  and  monitoring  of  breast  cancer 
therapy.  The  sensitivity  of  MRI/MRS  for  anatomical  delineation  is  very  high  and  the 
consensus  is  that  MRI  is  more  sensitive  in  detection  than  x-ray  mammography. 
Advantages  of  MRS  include  delivery  of  biochemical  information  about  tumor 
metabolism,  which  can  potentially  assist  in  the  staging  of  cancers  and  monitoring 
responses  to  treatment.  High  sensitivity  is  a  major  advantage  of  contrast  enhanced  MRI, 
but  its  diagnostic  relevance  in  the  future  will  largely  depend  on  improvements  in 
specificity.  Current  approaches  in  the  application  of  MRI  to  breast  tumors  aim  to 
improve  specificity  and  sensitivity  (4-17).  Increased  specificity  is  necessary  to  reduce  the 
number  of  biopsies  performed  to  confirm  false  positive  findings.  Diffusion-weighted 
imaging  (DWI)  is  another  MR  based  technique  that  probes  the  microstructure  of  tissues 
and  is  sensitive  to  the  degree  to  which  motion  of  water  molecules  is  restricted  in  relation 
to  how  packed  together  cells  are  (17,  18).  It  has  been  reported  that  high  resolution  DWI 
may  add  valuable  functional  information  to  conventional  MR  protocols  with  short 
measurement  times  for  the  diagnosis  of  breast  cancer  and  improve  the  specificity  of  MR 
imaging  (19-21).  However,  new  technological  developments  are  necessary  to  assess  their 
role  in  breast  diagnosis.  A  method  capable  of  identifying  biochemical  characteristics  non- 
invasively  in  the  tumor  lesions  that  can  be  used  in  conjunction  with  MRI  is  proton  ( 1 H) 
MR  Spectroscopy  (MRS).  Researchers  have  shown  that  'll  MRS  can  be  used  to 
characterize  breast  cancers  with  improved  diagnostic  accuracy  (22-26).  Unfortunately, 
multi-voxel  based  novel  MR  spectroscopic  imaging  (MRSI)  techniques  using  the  speed 
advantage  offered  by  echo-planar  imaging  (EPI)  and  improved  spectral  resolution  offered 
by  two-dimensional  (2D)  MR  spectroscopy  (MRS)  have  not  been  fully  explored  in  breast 
cancer  studies  so  far.  Hence,  a  major  task  of  this  project  is  to  combine  echo-planar 
correlated  spectroscopic  imaging  (EP-COSI)  with  DWI  approach  for  improving  the 
overall  specificity  of  breast  cancer  detection. 

Body: 

i)  Proposed  Task  1:  To  further  optimize  the  multi-voxel  based  extension  of  the  correlated 
spectroscopy  (COSY)  sequence,  in  which  two  spectral  encodings  will  be  combined  with 
two  spatial  encodings.  This  four-dimensional  (4D)  data  acquisition  scheme  will  be 
accomplished  utilizing  the  echo-planar  imaging  (EPI)  approach  that  is  commonly  used 
for  spatial  encoding  in  MRI  including  DWI.  (Months  1-6). 

This  task  was  completed  already  (September  2010-May  201 1)  as  reported  last  year. 

ii)  Proposed  Task  2:  To  evaluate  the  EP-COSI  data  using  a  breast  phantom  containing 
two  concentric  spheres,  the  inner  one  containing  several  metabolites  which  have  been 
reported  in  breast  tissues  surrounded  by  the  outer  phantom  containing  corn  oil  to  mimic 
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fatty  tissues  known  to  be  in  breast  tissues,  and  to  optimize  the  echo  speed-factor  and 
other  acquisition  parameters  using  the  phantom  (Months  1-6). 

This  task  was  completed  already  (September  2010-August  2011)  as  reported  again  last 
year. 


Hi)  Proposed  Task  3:  To  develop,  evaluate  and  optimize  the  prior-knowledge  basis  set 
spectra  using  the  GAMMA-simuIation  and  breast  phantom  solutions  as  prior  knowledge 
for  the  multi-voxel  based  COSY  spectra  recorded  using  the  3T  MRI  scanner  (Months  3- 

9f 

Using  the  GAMMA  library  (27),  a  prior-knowledge  basis-set  containing  metabolites  and 
lipids  was  prepared  (April  2010-August  201 1)  as  reported  last  year. 

iv)  Proposed  Task  4:  To  record  the  EP-COSI  spectra  in  the  fatty,  glandular  and  ductal 
areas  of  healthy  breasts.  Twenty  healthy  female  volunteers  (25-70  years  old)  with  no 
previous  history  of  breast  cancer  will  be  investigated.  (Months  6-24). 

Five  healthy  women  (age=36-59  years)  were  investigated  using  the  DWI  and  EP-COSI 
MRI  protocol  as  reported  earlier  (March  2011-July  2011).  Six  more  healthy  women 
(age=26-46  years)  without  any  prehistory  of  breast  cancer  have  been  studied  during  the 
2nd  year.  The  ADC  image  derived  from  the  DWI  data  recorded  in  a  34  yo  healthy  woman 
is  shown  in  Fig.  1.  Feasibility  of  recording  2D  COSY  spectra  in  multiple  regions  using 
the  EP-COSI  sequence  is  demonstrated  in  Fig.2.  Multi-voxel  metabolite  spectra  of  the 
fat  peak  at  5.3ppm  is  shown  in  Fig.2B  (28).  2D  MR  spectra  were  recorded  in  multiple 
regions  including  the  fatty  and  glandular  regions  of  the  34  yo  healthy  subject  as  shown  in 
Fig.2C  and  D.  C  and  D  were  extracted  COSY  spectra  from  the  fatty  and  glandular  regions 
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with  a  volume  of  1x1x2cm  .  These  findings  were  reproduced  in  other  10  ten  healthy 
women  as  shown  in  our  earlier  manuscript  using  a  single  voxel  (SV)  localized  L-COSY 
sequence  (29). 


Figure  1.  A)  An  axial  ADC  slice  image  recorded  in  a  34  yo  healthy  subject  using  the 
DWI  sequence. 
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Figure  2.  A)  An  axial  MRI  slice  image  recorded  in  the  same  healthy  subject  as  above 
showing  the  EP-COSI  localization;  the  small  white  box  shows  the  regions  of  interest 
(ROI)  localized  the  EP-COSI  sequence  and  the  large  white  box  with  16x16  grids  shows 
the  spatial  localization  encoded  by  the  EP-COSI  sequence.  B)  Multi-voxel  display  of  the 
EP-COSI  data  showing  the  olefenic  diagonal  peaks.  C)  and  D)  are  2ml  2D  COSY  spectra 
from  the  fatty  and  glandular  regions  extracted  from  the  EP-COSI  data. 

v)  Proposed  Task  5:  To  record  the  multi-voxel-based  2D  spectra  in  patients  with  benign 
and  malignant  breast  cancer.  The  breast  metabolite  and  lipid  concentrations  calculated 
from  the  multi-voxel  data  using  the  ProFit  algorithm  will  be  compared  with  LC-Model 
processed  ID  spectral  based  MRSI  data.  Twenty  patients  with  biopsy-proven  breast 
cancer  (ductal  carcinoma  and  invasive  lobular  cancer),  twenty  patients  with  benign 
breast  tumor  (fibroadenoma,  proliferative  fibrocystic  change  and  papillomas)  will  be 
investigated  (Months  6-24). 

Three  women  with  benign  (age=28-44  years)  and  2  with  malignant  breast  cancer 
(age=61  years)  have  been  investigated  using  the  4D  EP-COSI  sequence  on  the  3T  MRI 
scanner  using  a  dedicated  breast  MRI  coil.  Shown  in  Fig.  3  are  the  following:  A)  An  axial 
MRI  slice  image  recorded  in  a  44yo  benign  cancer  patient  showing  the  EP-COSI 
localization.  B)  Multi-voxel  EP-COSI  spectra  and  C)  a  2ml  COSY  spectrum  recorded 
from  the  affected  mass. 
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cancer  showing  the  EP-COSI  localization;  the  small  white  box  shows  the  regions  of 
interest  (ROI)  localized  the  EP-COSI  sequence  and  the  large  white  box  with  16x16  grids 
shows  the  spatial  localization  encoded  by  the  EP-COSI  sequence.  (B)  multi-voxel  2D 
COSY  spectra  covering  fatty  and  glandular  breast  regions.  (C)  A  2ml  benign  lesion- 
indicating  spectrum  showing  slightly  increased  water  and  slightly  reduced  lipids. 
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Figure  4.  (A)  An  axial  MRI  slice  image  recorded  in  a  61  yo  subject  with  malignant 
breast  cancer  (a  hypo-intense  mass)  showing  the  EP-COSI  localization;  the  small  white 
box  shows  the  regions  of  interest  (ROI)  localized  the  EP-COSI  sequence  and  the  large 
white  box  with  16x16  grids  shows  the  spatial  localization  encoded  by  the  EP-COSI 
sequence.  (B)  multi-voxel  2D  COSY  spectra  covering  fatty  and  glandular  breast  regions. 
(C)  and  (D)  show  1ml  spectra  showing  increased  water  and  slightly  reduced  lipids. 


Fig.  4  shows  the  following:  A)  An  axial  MRI  slice  image  recorded  in  a  61yo 
malignant  breast  cancer  patient  showing  the  EP-COSI  localization.  B)  Multi-voxel  EP- 
COSI  spectra;  C)  and  D)  show  2ml  COSY  spectra  recorded  from  the  affected  mass  and  a 
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healthy  region.  These  results  are  also  in  agreement  with  our  earlier  data  on  the  SV-based 
L-COSY  sequence  (29-30). 


vi)  Proposed  Task  6:  To  record  multi-slice  DWI  in  twenty  patients  with  biopsy-proven 
breast  cancer,  twenty  patients  with  benign  breast  tumor  and  tw>enty  healthy  women  and  to 
calculate  the  ADC  maps.  (Months  6-24). 

Figure  5.  An  axial  ADC  slice  image  recorded  in  the  44  yo  subject  with  a  benign  mass 
using  the  DWI  sequence. 


Figure  6.  An  axial  ADC  slice  image  recorded  in  the  61  yo  subject  with  a  malignant  mass 
in  the  left  breast  using  the  DWI  sequence. 
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The  diffusion  weighted  imaging  (DWI)  has  been  recorded  in  three  women  with  benign 
and  2  with  malignant  breast  cancer  3T  MRI  scanner  using  a  dedicated  breast  MRI  coil. 

vii)  Proposed  Task  7:  To  correlate  the  EP-COSI  findings  with  that  of  DWI  in 
differentiating  benign  from  malignant  breast  cancers,  and  to  calculate  specificity, 
sensitivity  and  accuracy  of  the  MRSI  and  DWI  data  in  differentiating  benign  from 
malignant  tumors.  (Months  6-24). 

This  task  will  be  completed  during  the  next  year. 

Key  Research  Accomplishments 

•  The  4D  EP-COSI  has  been  successfully  evaluated  in  1 1  healthy  women,  3  benign 
and  2  malignant  breast  cancer  patients  so  far  using  the  UCLA  Radiology  Siemens 
3T  MRI  scanner  equipped  with  a  dedicated  breast  phased-array  assembly.  As 
summarized  in  the  last  year  report,  this  sequence  is  available  now  at  UCLA  only 
and  is  not  supplied  by  any  of  MRI  manufacturers. 

•  The  DWI-MRI  protocol  was  successfully  evaluated  in  a  total  of  1 1  healthy 
women,  3  benign  and  2  malignant  breast  cancer  patients  so  far. 

•  We  have  demonstrated  that  the  EP-COSI  spectroscopic  imaging  sequence  can  be 
combined  with  a  clinical  breast  DWI  protocol  with  the  total  duration  of  less  than 
an  hour.  The  protocol  is  safe  enough  to  be  included  in  any  MRI  protocol  to  be 
evaluated  in  breast  cancer  for  improving  the  overall  specificity. 

•  We  are  currently  testing  retrospective  Maximum  Entropy  and  Compressed 
Sensing  of  the  4D  EP-COSI  data  so  that  the  acquisition  can  be  accelerated  in  the 
near  future  to  5-10  minutes  which  will  lead  to  less  inconvenience  to  patients 
suffering  from  breast  cancer. 


Reportable  Outcomes: 

A.  Peer-reviewed  Publications:  Burns  B,  Wilson  N  and  Thomas  MA.  Maximum 

Entropy  Reconstruction  of  Non-Uniformly  Under-Sampled  Multi-dimensional  MR 

Spectroscopic  Imaging  in  vivo.  Magn  Reson  Med  2012  (submitted). 

B.  Presentations:  There  were  three  conference  presentations  containing  four  abstracts. 

•  A  Poster  presentation  at  the  International  Society  of  Magnetic  Resonance  in  Medicine 
(ISMRM)  Conference  2012,  Melbourne,  Australia.  Maximum  Entropy  Reconstruction 
of  Correlated  Spectroscopy  of  Human  Breast  in  vivo.  May  2012 

•  A  Poster  presentation  at  the  International  Society  of  Magnetic  Resonance  in  Medicine 
(ISMRM)  Conference  2012,  Melbourne,  Australia.  Maximum  Entropy  Based 
Reconstruction  of  Echo-Planar  Correlated  Spectroscopic  Imaging  of  Human  Breast 
in  vivo.  May  2012 

•  A  Poster  presentation  at  the  International  Society  of  Biomedical  Engineering  (ISBI) 
Conference  2012,  Barcelona,  Spain.  Maximum  Entropy  Reconstruction  of  Echo- 
Planar  Correlated  Spectroscopic  Imaging  of  Human  Breast  In  Vivo.  May  2012. 
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•  A  Poster  presentation  at  the  5 1st  Experimental  NMR  Conference  (ENC)  ,  Orlando, 
Florida.  Accelerated  Multi-voxel  two-dimensional  in  vivo  spectroscopy  using 
Compressed  Sensing.  April  2012 


C.  Books:  None  on  Breast  Cancer  Research  based. 

Conclusions:  The  scanning  protocol  including  EP-COSI  and  DWI-MRI  has  been 
evaluated  in  1 1  healthy  women,  three  benign  and  2  malignant  breast  cancer  so  far.  We 
will  continue  to  recruit  10  malignant,  10  benign  and  5  healthy  women  during  the  next 
year.  Our  inability  in  recruiting  more  breast  cancer  patients  has  been  discussed  with  the 
UCLA  mammo  co-investigators,  namely  Dr.  Naneltte  DeBruhl  and  Dr.  Lawrence 
Bassett.  They  have  promised  us  to  improve  the  patient  recruitment  in  the  3rd  year. 
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Introduction 

4D  Echo-Planar  based  Correlated  Spectroscopic  Imaging  (EPCOSI) 
allows  for  the  simultaneous  acquisition  of  two  spatial  (ky,  kx)  and  two 
spectral  (t2,  tx)  dimensions  in  vivo  in  a  single  Magnetic  Resonance  (MR) 
experiment  [1]. 

The  EPI  readout  accelerates  the  direct  acquisition  of  the  kx  and  t2 
dimensions  but  the  ky  and  t2  dimensions  are  incrementally  acquired  as 
indirect  dimensions. 

Non-uniform  under-sampling  (NUS)  the  ky  and  t2  dimensions  to  25% 
and  reconstructing  the  missing  data  with  Maximum  Entropy  (MaxEnt) 
preserves  the  spectral  and  spatial  resolution  of  the  acquired  data  and 
reduces  the  scan  time  from  20  to  5  minutes  in  healthy  human  breast  in 
vivo. 

Materials  and  Methods 

An  EPCOSI  scan  of  healthy  breast  with  the  following  parameters: 
1x1x2cm3  voxel  size,  TR/TE=1.5s/30ms,  50  t2  increments,  16x16cm2  FOV 
was  retrospectively  NUS  to  25%  of  the  ky-^  plane  using  the  mask  in  Fig. 
2. 

Missing  data  in  the  ky-t2  plane  was  reconstructed  using  MaxEnt  [2,3]: 
maximize  S1/2(f)  s.t.  C(f)  <  =  a  (1) 
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Figure  1  (Top)  Selected  2D  COSY  spectra  from  a  4D  EPCOSI  scan  of  healthy,  fatty  breast  highlighted  in  red  in  the  MRI 
localizer  image  for  Al)  fully  sampled  Bl)  25%  NUS  Cl)  25%  NUS  with  MaxEnt  reconstruction.  (Bottom)  Spatial 
distribution  of  olefinic  cross-peaks  of  unsaturated  fatty  acids  highlighted  in  light  blue  for  A2)  fully  sampled  B2)  25% 
NUS  C2)  25%  NUS  with  MaxEnt  reconstruction.  _ 

Algorithm  :  Cambridge  MaxEnt  (NUS  data,  sample  density) 


where  f  is  the  estimated  fully-sampled  spatially  distributed  spectra,  C(f)  is  the  time-domain  fidelity  constraint,  a  is  the 
noise  standard  deviation,  and  S1/2(f)  is  the  spin-1/;  entropy  of  the  estimated  spectra  [4]. 

The  Cambridge  MaxEnt  algorithm  uses  a  variant  of  the  conjugate  gradient  method  to  solve  (1)  as  shown  in  Fig.  3. 

Results  and  Discussion 

Al  is  a  voxel  from  a  fully  sampled  4D  EPCOSI  scan  of  healthy,  fatty  breast  in  vivo.  Bl  is  the  same  voxel  after  25%  NUS  of 
the  ky-t2  plane  and  has  significant  artifacts  from  the  large  diagonal  peak  at  (2.2ppm,  2.0ppm)  obscuring  the  much  smaller 
cross-peak  highlighted  in  blue.  Cl  shows  the  MaxEnt  reconstruction  of  Bl  and  has  equivalent  spectral  peaks  to  the  fully 
sampled  data  in  Al  with  the  artifacts  removed. 

A2,  B2,  and  C2  show  the  spatial  distribution  of  the  olefinic  cross  peaks  highlighted  in  blue  for  the  fully  sampled,  25% 

NUS,  and  25%  NUS  with  MaxEnt  reconstruction  from  a  4D  EPCOSI  scan  of  healthy,  fatty  breast.  B2  shows  significant 
spatial  incoherent  artifacts  from  the  NUS  that  are  removed  by  MaxEnt  reconstruction  in  C2. 

The  normalized  Root  Mean  Square  Error  (nRMSE)  of  B1/B2  is  .0047  and  .0014  for  C1/C2. 

Conclusion 

We  have  shown  it  is  possible  to  simultaneously  NUS  a  spectral  and  a  spatial 
dimension  of  an  in  vivo  EPCOSI  scan  to  25%  and  successfully  reconstruct  the  data. 
This  acceleration  translates  into  a  clinically  viable  5  minute  EPCOSI  scan  time. 
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Figure  2  Poisson  distributed  25%  sample 
mask  for  ky  and  tx  dimensions.  Sample 
density  follows  sin-modulated  tx  and 


High  level  overview  of  MaxEnt: 
repeat 

Calculate  gradients  and  hessians: 
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Create  local  -,“l  order  approximation  from  basis: 


(Hi) 


I  SO,)  -  S(f)  +  y’b’g,  -  Wb’H,  by 
|  C(y)  =  C(f)  +  y’b’g,  +  —y'b’U,by 


Rotate  basis  to  make  subspace  Cartesian; 

Iy=>X,  b'H ,  b=>  /,  bmUc  b  =>  Mdiag 
S(x)  =  S/  +  x’b’g,  -  ti'/i 

C(r)  -  C,  +  x’b'g,  +  xx’Ma,gx 
Find  x  that  minimizes  C(x)  and  maximizes  S(x)  from  (iv): 
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Introduction  -  The  altered  metabolism  of  breast  cancer  tissues  gives  rise  to  changes  in  metabolite  concentrations  that  are  detectable  using  two- 
dimensional  (2D)  Magnetic  Resonance  Spectroscopy  (MRS)  [1].  2D  Localized-Correlated  Spectroscopy  (L-COSY)  MRS  has  been  shown  to  increase 
the  specificity  and  sensitivity  of  tumor  grade  classification  when  used  with  traditional  techniques,  such  as  Dynamic  Contrast  Enhanced  (DCE)-MRI 
[1],  However,  because  of  the  number  of  ti  increments  used  to  construct  the  indirect  spectral  dimension  (F  i),  the  L-COSY  sequence  requires  on  the 
order  of  25  minutes  for  100  q  increments  and  8  averages  [2],  To  reduce  scan  times  to  a  clinically  acceptable  12  minutes,  only  45  q  increments  are 
typically  used,  but  this  reduces  spectral  resolution  along  F  | .  Maximum  entropy  (MaxEnt)  image  reconstruction  techniques  have  been  used  in  Nuclear 
Magnetic  Resonance  (NMR)  to  non-unifonnly  under-sample  (NUS)  indirect  spectral  dimensions  and  then  reconstruct  the  most  statistically  likely 
fully-sampled  multi-dimensional  spectrum  [3,4],  This  technique  can  be  used  to  accelerate  the  collection  of  2D  L-COSY  data  in  vivo  by  selectively 
under-sampling  along  q.  We  show  that  under-sampling  the  indirect  dimension  by  ~50%  and  reconstructing  the  spectrum  using  MaxEnt  preserves  the 
metabolite  diagonal  peaks:  olefmic  fat  (UFD),  methyl  fat  (FMETD),  and  choline  (Cho)  and  the  cross  peaks:  unsaturated  fatty  acid,  right  (UFR), 
unsaturated  fatty  acid,  left  (UFL),  and  triglyceryl  fat  (TGFR)  which  have  demonstrated  clinical  importance. 

Methods  -  To  qualitatively  measure  the  performance  of  the  MaxEnt  reconstruction  algorithm  on  NUS  in  vivo  data,  fully  sampled  2D  L-COSY  scans 
of  both  healthy  (n=10)  and  malignant  (n=8)  breast  tissues  were  acquired  on  a  3T  Siemens  Trio  scanner  using  the  following  parameters:  lcm3  voxel 
size,  TR/TE=2s/30ms,  45  q  increments,  8  averages,  and  no  water  suppression.  Both  scans  were  retrospectively  under-sampled  to  22  lines  along  the  q 
dimension,  apodized  along  the  q  and  t2  dimensions,  and  then  reconstructed  using  the  MaxEnt  algorithm  to  45  lines.  MaxEnt  is  a  constrained  convex 
optimization  algorithm  that  uses  a  variant  of  the  conjugate  gradient  method  to  iteratively  solve  the  inverse  problem  [5]: 

maximize  S i/2(f)  s-t-  -  D\\2  <=  a  (1) 

where /is  the  estimated  fully-sampled  spectrum  at  each  iteration,  F1  is  the  inverse  Fourier  transform,  K  is  the  NUS  matrix,  D  is  the  time-domain 
measured  data,  0  is  the  noise  standard  deviation,  and  S  112(f)  is  the  spin-!4  entropy  of  the  estimated  spectrum  [3],  The  spectrum  with  the  highest 
entropy  is  that  which  conforms  to  the  uniform  distribution  and  has  a  flat  baseline.  Therefore,  by  maximizing  the  spin-VS  entropy  in  (1),  MaxEnt 
enforces  sparsity  of  the  estimated  spectrum  in  the  frequency  domain  and  reduces  incoherent  aliasing  artifacts  in  the  spectrum  caused  by  NUS.  The 
fidelity  constraint  ensures  that  peaks  in  the  estimated  spectrum  must  come  from  the  sampled  data  to  within  a  tolerance  of  the  noise.  The 
reconstruction  is  performed  over  both  dimensions  simultaneously  as  opposed  to  a  series  of  1 D  reconstructions  as  implemented  elsewhere  [4], 


Results  and  Discussion  -  Figures  1 A  and  IB  show  the  2D  L-COSY  spectra  of  the  fully  sampled  and  ~50%  MaxEnt  reconstruction  of  the  NUS 
healthy  breast  tissue,  respectively.  As  can  be  seen,  the  major  diagonal  and  cross  peaks  present  in  healthy  fatty  breast  tissue  are  present  in  both 
spectra.  The  peak  locations,  line  shapes,  and  amplitudes  of  the  NUS  data  are  roughly  equivalent  to  the  fully  sampled  data;  however  a  slight  loss  in 
spectral  resolution  along  F ,  is  evident  by  the  elongation  of  the  UFR  peak  in  the  MaxEnt  reconstruction.  There  is  no  corresponding  loss  in  resolution 
along  F2.  Figures  2A  and  2B  show  2D  L-COSY  spectra  of  the  fully  sampled  and  the  ~50%  reconstruction  of  the  NUS  malignant  fatty  breast  tissue, 
respectively.  The  presence  of  water  in  fatty  breast  tissue  is  an  indicator  of  breast  cancer,  and  because  of  its  high  concentration,  its  broad  line-width 

obscures  many  of  the  cross-peaks  seen  in  Figure  1.  However, 
the  UFR  cross-peak  is  clearly  visible  in  figures  2A  and  2B,  and 
the  MaxEnt  reconstructed  peak  is  comparable  to  the  fully 
sampled  peak,  even  in  the  presence  of  noise.  The  Cho  diagonal 
peak  in  figure  2,  which  is  a  sign  of  increased  metabolic  activity 
and  malignancy,  is  fully  resolved  in  the  MaxEnt  reconstruction 
and  shows  good  agreement  with  the  peak  in  the  fully  sampled 
spectrum. 


Conclusions  -  We  have  shown  that  it  is  possible  to  under¬ 
sample  the  L-COSY  sequence  by  up  to  ~50%  and  reconstruct  in 
vivo  spectra  from  as  few  as  22  q  lines  that  show  similar  spectral 
characteristics  to  spectra  from  fully  sampled  data.  This 
acceleration  translates  into  a  sub-6  minute  scan  time  if  8 
averages  are  used  and  increases  the  potential  use  of  MRS  for 
breast  cancer  screening  in  the  clinic.  If  increased  spectral 
resolution  is  needed,  the  MaxEnt  reconstruction  is  not  limited 
to  45  q  points  and  can  reconstruct  many  more  NUS  points  than 
what  was  used  in  these  experiments. 


Fig  2:  in  vivo  2D  L-COSY  spectra  of  a  25 
year  old  healthy  fatty  breast  tissue  (A)  fully 
sampled  along  ti  (B)  50%  reconstruction 
along  ti 


Fig  2:  in  vivo  2D  L-COSY  spectra  of  a  55 
year  old  malignant  tumor  in  fatty  breast  tissue 
(A)  fully  sampled  along  t,  (B)  50% 
reconstruction  along  ti 
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Introduction  -  2D  Localized  Correlated  Spectroscopy  (L-COSY)  has  been  shown  to  be  a  powerful  tool  in  the  detection  of  breast  cancer  but  is 
limited  to  a  single  voxel  [1],  The  recently  introduced  Echo-Planar  based  Correlated  Spectroscopic  Imaging  (EPCOSI)  sequence  allows  for  the 
simultaneous  acquisition  of  two  spatial  (ky,  kx)  and  two  spectral  (t2,  /)  dimensions  in  a  single  experiment  [2],  The  4D  EPCOSI  sequence  interleaves 
the  acquisition  of  the  kx  and  t2  dimensions  along  the  EPI  read-out,  while  ky  and  /  lines  are  incrementally  acquired  as  indirect  dimensions,  requiring 
20  minutes  for  a  typical  scan.  To  reduce  scan  times  using  standard  Fast  Fourier  Transform  (FFT)  based  reconstruction  techniques  would  require  the 
reduction  of  either  the  ky  spatial  or  tj  spectral  dimensions,  and  a  corresponding  reduction  in  resolution.  Maximum  entropy  (MaxEnt)  image 
reconstruction  techniques  have  been  used  in  Nuclear  Magnetic  Resonance  (NMR)  to  non-uniformly  under-sample  (NUS)  indirect  spectral 
dimensions  and  in  imaging  to  NUS  spatial  dimensions  then  reconstruct  the  fully-sampled  multi-dimensional  spectra  or  image  [3,4,5],  This  technique 
can  be  used  to  accelerate  the  collection  of  4D  EPCOSI  data  in  vivo  by  selectively  under-sampling  along  ky  and  /.  We  show  that  under-sampling  the 
indirect  dimensions  by  ~25%  and  reconstructing  the  spatially  localized  EPCOSI  spectra  using  MaxEnt  preserves  the  spatial  distribution  of  the 
metabolite  diagonal  peaks  found  in  healthy  breast  tissue  and  reduces  the  scan  time  to  5  minutes. 

Methods  -  An  EPCOSI  scan  of  healthy  breast  was  perfonned  on  a  3T  Siemens  Trio  scanner  with  the  following  parameters: 

1x1x2cm3  voxel  size,  TR/TE=1.5s/30ms,  50  /  increments,  16xl6cnr  FOV,  1  average  with  water  suppression  and  2  averages 
without  water  suppression.  The  scan  was  retrospectively  under-sampled  to  ~25%  of  the  ky-ti  plane,  using  the  map  in  Fig.  1. 

The  NUS  spatial  and  spectral  dimensions  were  simultaneously  reconstructed  using  MaxEnt.  The  MaxEnt  algorithm  uses  a 
variant  of  the  conjugate  gradient  method  to  iteratively  solve  the  constrained  convex  optimization  problem  [5]: 

maximize  Si /z(f)  s.t.  \\RlKf- D\\2  <=  a  (1) 

where /is  the  estimated  fully-sampled  spatially  distributed  spectra  at  each  iteration,  R1  is  the  inverse  Fourier  transform,  K  is 
the  NUS  matrix,  D  is  the  time-domain  measured  data,  a  is  the  noise  standard  deviation,  and  Si^ff)  is  the  spin- 14  entropy  of  the 
estimated  spectra  [3].  The  spatially  distributed  spectra  with  the  highest  entropy  are  those  which  conform  to  the  unifonn 
distribution  and  have  a  flat  baseline.  Therefore,  by  maximizing  the  spin- 14  entropy  in  (1),  MaxEnt  enforces  sparsity  of  the 
estimated  spectra  in  the  frequency  and  k-space  domains  and  reduces  incoherent  aliasing  artifacts  in  the  spectra  caused  by 
NUS.  The  fidelity  constraint  ensures  that  peaks  in  the  estimated  spectra  and  their  spatial  distribution  must  come  from  the 
sampled  data  to  within  a  tolerance  of  the  noise.  The  reconstruction  is  performed  over  all  dimensions  simultaneously  as 
opposed  to  a  series  of  ID  reconstructions  as  implemented  elsewhere  [5], 

Results  and  Discussion  -  Figures  2  A  and  2B  show  the  fully  sampled  and  25%  MaxEnt  reconstruction  of  the  selected  2D 
COSY  spectra  from  a  voxel  of  healthy  fatty  breast  tissue.  As  can  be  seen,  the  spectral  characteristics  of  the  NUS  data  with 
MaxEnt  reconstruction  are  equivalent  to  the  fully  sampled  data;  the  same  diagonal  and  cross  peaks  are  present  in  both  spectra 
without  a  significant  increase  in  noise.  Figures  3A  and  3B  show  the  spatial  distribution  of  the  olefinic  cross  peak  selected  in 
Figure  2  for  the  fully  sampled  and  25%  MaxEnt  reconstructed  4D  EPCOSI  of  healthy  fatty  breast  tissue.  The  spatial  distributions  between  the  fully 

sampled  and  MaxEnt  reconstruction  show  very  good  agreement, 
with  the  exception  of  some  spatial  bleeding  on  the  outer  edges 
of  the  distribution. 

Conclusions  -  We  have  shown  that  it  is  possible  to  under¬ 
sample  both  the  spectral  and  spatial  dimensions  of  an  in  vivo 
EPCOSI  sequence  by  up  to  25%  and  reconstruct  the  spectra 
with  similar  spatial  distributions  and  spectral  characteristics  to 
the  fully  sampled  data.  This  acceleration  translates  into  a  5 
minute  EPCOSI  scan  time  and  increases  the  potential  use  of  4D 
MRS  in  the  clinic.  Compressed  Sensing  (CS)-based 
reconstruction  can  be  use  alternatively  for  processing  the  NUS 
EPCOSI  data.  Further  research  regarding  the  merits  and 
demerits  of  the  two  reconstruction  methods  is  in  progress. 
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Figure  1:  4x 
under-sampling 
mask  for  ky  and  ti 
dimensions.  White 
points  are  sampled. 


Figure  2:  Selected  2D  COSY  spectra  of 
healthy  fatty  breast  tissue  highlighted  in 
red  for  A)  fully  sampled  and  B)  4x  under¬ 
sampled  with  MaxEnt  reconstruction. 


Figure  3:  Spatial  distribution  of  Olefinic 
cross-peaks  highlighted  in  purple  in  Fig.  2 
for  A)  fully  sampled  and  B)  4x  under¬ 
sampled  with  MaxEnt  reconstruction. 
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Abstract 

Purpose:  To  reduce  scan  times  of  4D  Magnetic  Resonance  Spectroscopic  Imaging  by 
applying  non-uniform  under-sampling  (NUS)  to  the  indirectly  acquired  dimensions  (ky 
and  ti)  and  reconstructing  those  dimensions  using  Maximum  Entropy  (MaxEnt)  iterative 
reconstruction. 

Methods:  Simulated  and  in  vivo  4D  Echo-Planar  based  Correlated  Spectroscopic 
Imaging  (EP-COSI)  data  were  retrospectively  under-sampled,  and  in  vivo  4D  Echo- 
Planar  J-Resolved  Spectroscopic  Imaging  (EP-JRESI)  prostate  data  was  prospectively 
undersampled  in  the  ky-ti  plane.  The  under-sampled  data  sets  were  reconstructed  using 
the  Cambridge  iterative  algorithm  to  solve  the  MaxEnt  reconstruction  problem. 

Results:  Simulated  data  reconstructions  demonstrated  the  MaxEnt  reconstruction 
improves  data  fidelity  considerably  over  the  under-sampled  data  at  SNRs  from  2  to  20 
and  percent  under-samplings  of  20  to  80%.  25%  NUS  was  applied  retrospectively  to  EP- 
COSI  and  prospectively  to  EP-JRESI  in  vivo  data.  The  NUS  EP-COSI  and  EP-JRESI  in 
vivo  data  were  under-sampled  to  25%  and  the  incoherent  artifacts  caused  by  the  under¬ 
sampling  were  successfully  removed  by  the  MaxEnt  reconstruction. 

Conclusion:  MaxEnt  reconstruction  preserved  the  spectral-spatial  resolution  of  the  non- 
uniformly  under-sampled  data  and  reduced  EP-JRESI  scan  times  from  24  to  6  minutes  in 
human  prostate  and  can  potentially  reduce  EP-COSI  scan  times  from  20  to  5  minutes  in 
human  breast. 

Keywords:  EP-COSI,  EP-JRESI,  Maximum  Entropy,  Non-Uniform  Under-Sampling, 
Spectroscopy,  Spectroscopic  Imaging 
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INTRODUCTION 

The  altered  metabolism  of  cancer  gives  rise  to  changes  in  metabolite  concentrations  that 
can  potentially  be  detected  non-invasively  using  one-dimensional  (ID)  magnetic 
resonance  spectroscopy  (MRS)  in  vivo  [1,2, 3, 4, 5].  However,  the  overlap  of  spectral 
peaks  in  ID  MRS  is  a  major  concern  for  identifying  individual  metabolites.  Two 
dimensional  (2D)  MRS  has  increased  spectral  dispersion  over  ID  MRS  and  can 
disentangle  over-lapping  complex  spectral  peaks  [6],  Single  voxel  2D  MRS  has  been 
shown  to  increase  the  specificity  and  sensitivity  of  tumor  grade  classification  when  used 
with  traditional  techniques,  such  as  dynamic  contrast  enhanced  (DCE)-MRI  in  the  breast 
[7]  and  the  Gleason  grading  system  in  the  prostate  [8].  However,  requiring  multiple  fi 
increments  per  voxel  to  fonn  the  second  spectral  dimension  limits  its  ability  to  provide 
multi-voxel  coverage  due  to  long  scan  times  needed  to  combine  two  spectral  and  two 
spatial  dimensions. 

With  the  advent  of  Echo-Planar  Spectroscopic  Imaging  (EPSI),  magnetic  resonance 
spectroscopic  imaging  (MRSI)  scans  with  one  spectral  dimension  could  be  completed 
within  clinically  acceptable  times  by  interleaving  the  acquisition  of  a  spatial  and  spectral 
dimension  within  the  EPSI  readout  [9]  [10]  [1  lj.The  recently  introduced  Echo-Planar 
based  Correlated  Spectroscopic  Imaging  (EP-COSI)  [12]  and  Echo-Planar  J-Resolved 
Spectroscopic  Imaging  (EP-JRESI)  [13]  sequences  allow  for  the  simultaneous  acquisition 
of  two  spatial  (ky,  kx)  and  two  spectral  (t2,  ti)  dimensions  in  a  single  recording  to  fonn  4D 
MRSI.  Both  sequences  interleave  the  acquisition  of  the  kx  and  t2  dimensions  within  the 
EPSI  read-out,  but  ky  and  fi  are  incrementally  acquired  as  indirect  dimensions  during 
each  TR.  The  EP-COSI  and  EP-JRESI  sequences  have  the  benefits  of  increased  spectral 
dispersion  and  multi-voxel  support;  however  their  scan  times  are  directly  proportional  to 
the  number  of  increments  in  the  ky  and  fi  dimensions  and  can  take  20  to  40  minutes  using 
typical  parameters  which  is  too  long  to  be  used  for  a  routine  clinical  protocol  [12]  [13]. 

When  using  the  Fast  Fourier  Transform  (FFT),  reducing  4D  EP-COSI  and  EP-JRESI 
scan  times  requires  the  reduction  of  either  the  ky  spatial  or  fi  spectral  dimensions  through 
truncation  or  lower  sampling  rates,  and  a  corresponding  unwanted  reduction  in  resolution 
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or  bandwidth  with  the  potential  for  aliasing.  However,  non-unifonn  under-sampling  of 
the  spatial-spectral  ky-ti  plane  in  combination  with  iterative  non-linear  reconstruction 
techniques  can  be  used  to  accelerate  the  collection  of  4D  MRSI  data  in  vivo  while 
preserving  the  spatial  and  spectral  resolution  [14]. 

Earlier  work  has  demonstrated  the  feasibility  of  under-sampling  the  mixed-domain  ky-ti 
plane  of  a  4D  EP-JRESI  data  set  and  reconstructing  the  missing  points  with  Compressed 
Sensing  (CS)  [14],  a  popular  method  of  under-sampled  iterative  image  reconstruction  that 
promotes  data  sparsity  in  the  reconstruction  domain  and  data  fidelity  in  the  sample 
domain  [15]  [16].  The  nature  of  spatial-spectral  incoherent  artifacts  in  the  ky-fi  plane 
were  explored,  and  it  was  shown  that  / /-norm- based  CS  reconstruction  is  a  viable  means 
of  reducing  the  scan  times  of  4D  EP-JRESI  in  vivo  through  under-sampling.  In  recent 
years,  CS  reconstruction  has  been  successfully  applied  to  under-sampled  MRI  [17]  [18], 
3D  MRSI  [19],  dynamic  MRI  [20]  [21],  and  a  number  of  other  promising  applications  not 
necessarily  related  to  medical  imaging  [22]  [23]. 

Maximum  Entropy  (MaxEnt)  image  reconstruction  is  an  alternative  non-linear  iterative 
reconstruction  technique  to  CS.  It  promotes  the  entropy  of  the  data  in  the  reconstruction 
domain  while  preserving  data  fidelity  in  the  sample  domain  [24,25].  MaxEnt  has  been 
successfully  applied  to  under-sampled  images  in  astronomy  and  multi-dimensional 
spectra  in  NMR  [25]  [26]  [27].  However,  MaxEnt  has  not  been  applied  to  the  mixed- 
domain  ky-ti  plane  of  a  4D  MRSI  data  set  in  vivo. 

In  information  theory,  entropy  is  a  measure  of  the  uncertainty  of  a  variable  state  and  the 
probability  of  its  occurrence;  the  more  uncertain  its  state,  the  higher  its  entropy.  States 
that  are  all  equally  likely  and  whose  probabilities  conform  to  the  unifonn  distribution 
have  the  maximum  possible  entropy  [28],  Entropy  in  MRSI  represents  a  measure  of  the 
phase  coherence  of  an  ensemble  of  spins  within  a  given  voxel  volume  [26].  Nonnalized 
for  the  concentration  of  spins  within  the  volume,  high-amplitude  NMR  signals  represent 
states  of  full  phase  coherence  and  low  entropy,  while  low  amplitude  NMR  signals 
represent  states  of  low  phase  coherence  and  high  entropy.  If  the  phase  distribution 
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conforms  to  a  state  of  low  phase  coherence,  such  as  a  uniform  distribution,  the  entropy  is 
maximized  as  it  is  in  information  theory.  Therefore,  by  maximizing  the  entropy  of  the 
spatial,  spectral-domain  in  MaxEnt,  under-sampling  artifacts  are  removed  from  the 
reconstruction  because  they  represent  states  of  high  phase  coherence  and  low  entropy  that 
are  not  present  in  the  sampled  data. 

In  this  paper,  we  show  that  MaxEnt  reconstruction  is  a  viable  technique  to  reduce  scan 
times  by  successfully  reconstructing  under-sampled  4D  MRSI  data.  We  quantitatively 
characterize  the  MaxEnt  reconstruction  by  showing  results  for  retrospectively  under¬ 
sampled  simulated  4D  EP-COSI  data  at  varying  levels  of  SNR  and  percent  under¬ 
sampling.  Additionally,  we  show  that  under-sampling  the  ky-ti  plane  to  25%  and  using 
MaxEnt  to  reconstruct  an  in  vivo  4D  EP-COSI  data  set  retrospectively  and  an  in  vivo  EP- 
JRESI  data  set  prospectively  preserves  the  spatial-spectral  resolution  of  the  spectra  and 
reduces  scan  times  to  5  minutes  in  the  breast  and  6  minutes  in  the  prostate  using  typical 
parameters. 

METHODS 
MR  Simulations 

The  effects  of  SNR  and  percent  under-sampling  on  the  MaxEnt  reconstruction  were 
quantitatively  assessed  using  a  noise-free  simulated  4D  EP-COSI  data  set  that  contained 
either  N-acetylaspartate  (NAA),  creatine  (Cr),  choline  (Cho),  lactate  (Lac),  or  nothing  in 
each  voxel,  as  represented  in  the  top  of  Figure  1  by  a  diagonal  peak  from  each  metabolite. 
Each  metabolite  was  simulated  using  the  GAMMA  NMR  libraries  [29]  from  a  Localized 
2D  Correlated  Spectroscopy  (L-COSY)  sequence  [30]  with  the  following  parameters:  100 
ti  increments,  1024  points  in  t2,  TR/TE=  1.5s/30ms,  and  spectral  bandwidths  of  1250Hz 
and  2000Hz  along  Fi  and  F2,  respectively.  Each  2D  spectrum  was  line -broadened  by  1  Hz 
and  apodized  by  a  sine-squared  filter  along  ti  and  a  skewed  sine-squared  filter  with  skew 
parameter  0.5  along  t2.  They  were  then  copied  into  an  8x8  spatial  grid  to  simulate 
spatially  distributed  metabolites  as  follows:  the  upper  left  quadrant  contained  2x2  voxels 
of  NAA,  the  upper  right  quadrant  contained  2x2  voxels  of  Cr,  the  lower  left  quadrant 
contained  2x2  voxels  of  Cho,  and  the  lower  right  quadrant  contained  2x2  voxels  of  Lac. 
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As  a  result  of  under-sampling  the  ky-tj  plane,  the  spatial-spectral  artifacts  caused  the 
NAA  and  Cho  voxels  to  alias  into  each  other  and  the  Cr  and  Lac  voxels  to  alias  into  each 
other,  as  shown  in  the  bottom  of  Figure  1.  This  changed  the  NAA  to  Cho  and  Lac  to  Cr 
ratios  as  the  spatial  and  spectral  separation  between  the  metabolites  broke  down. 

The  simulated  4D  EP-COSI  data  set  was  retrospectively  reconstructed  using  MaxEnt  at 
20,  40,  60  and  80%  under-sampling  and  SNRs  of  2  to  20  in  increments  of  2.  The  SNR 
was  varied  by  simulating  different  levels  of  thermal  noise  in  the  data  set  by  adding 
Gaussian  noise  in  quadrature  to  the  noise-free  4D  EP-COSI  data  set  as  in  [3 1].  The 
desired  SNR  was  achieved  by  ensuring  the  Root  Mean  Square  (RMS)  of  the  additive 
noise  was  equal  to  1/SNR  of  the  noise-free  data  set  RMS,  such  that: 

RMS(Noise  Free ) 

noisy  data  =  noise  free  -\ - 7 - - - *  Noise  fl] 

7  J  RMS  (Noise)  *  SNR  L  J 

Because  the  additive  noise  was  random,  each  SNR  was  simulated  and  reconstructed  20 
times  per  sample  mask  to  account  for  random  fluctuations  in  the  reconstruction.  The 
sampling  masks  were  created  using  the  2D  Poisson-gap  method  described  below. 

MR  Spectroscopic  Imaging 

A  4D  EP-COSI  scan  of  a  34  year  old  healthy  human  breast  was  acquired  on  a  Siemens 
3T  Trio  scanner  with  the  following  parameters:  lxlxl  cm'  voxel  size,  50  ti  increments, 
TR/TE/averages  =  1.5s/30ms/l,  a  16x16cm2  FOV,  and  spectral  bandwidths  of  1250Hz 
and  1 190Hz  along  Fi  and  F2,  respectively.  The  scan  was  fully  sampled  and  took  20 
minutes  to  complete.  Prior  to  being  under-sampled  and  reconstructed,  coil  combination 
and  eddy  current  correction  were  applied  to  the  fully  sampled  data  set.  It  was  then 
apodized  by  a  sine-squared  fdter  along  ti  and  a  skewed  sine-squared  fdter  with  skew 
parameter  0.5  along  t2.  The  fully  sampled  data  set  was  then  retrospectively  under¬ 
sampled  to  25%  in  the  ky-ti  plane  using  the  mask  shown  in  Figure  4  and  reconstructed  by 
MaxEnt. 
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The  prospective  MaxEnt  reconstruction  was  performed  on  a  4D  EP-JRESI  scan  of  a  7 1 
year  old  human  prostate  with  cancerous  tumors  in  the  left  and  right  base  (Gleason  score 
of  3+3,  prostatic  specific  antigen  of  8).  The  scan  was  acquired  on  a  Siemens  3T  Trio 
scanner  using  the  single-channel  endorectal  coil  with  the  same  parameters  as  the  EP- 
COSI  breast  scan  except  64  fi  increments  and  a  spectral  bandwidth  of  ±500Hz  along  Fi 
was  used.  The  ky-fi  mask  shown  in  Figure  5  was  used  during  the  scan  to  acquire  25%  of 
the  ky-ti  plane  during  the  acquisition.  Eddy  current  correction  was  applied  to  the  4D 
prostate  data  after  MaxEnt  reconstruction. 

Sample  Mask  Generation 

The  2D  Poisson-gap  sample  masks  used  were  derived  from  ID  Poisson-gap  sample 
masks  [32].  Poisson-gap  sample  masks  avoid  large  gaps  between  samples,  which  are 
detrimental  to  the  reconstruction,  while  ensuring  the  samples  are  randomly  distributed 
[33]. 

The  effects  of  the  sample  mask  on  the  SNR  and  line-shape  of  spectral  reconstructions  are 
well  documented  [33]  [34]  [35].  Sample  masks  that  follow  the  time-domain  NMR  signal 
envelope  and  sample  more  points  at  higher  SNR  have  lower  RMS  errors  (RMSE)  and 
non-linearity  compared  to  sample  masks  that  do  not.  Because  the  fi  dimension  in  a 
filtered  EP-COSI  scan  has  a  sine-bell  signal  envelope  and  the  fi  dimension  in  an  EP- 
JRESI  scan  has  an  exponentially  decaying  signal  envelope,  the  2D  Poisson-gap  sample 
masks  were  modulated  along  fi  with  sine  or  exponential  decay  functions,  respectively. 
The  ky  dimension  for  both  scan  types  was  modulated  by  an  exponential  decay  function 
similarly  to  what  has  been  used  previously  to  maximize  spatial  SNR  [36]. 

The  rate  parameter,  A,  determines  both  the  mean  and  variance  of  a  Poisson  distribution. 
The  probability  of  generating  a  sample  gap,  g,  from  a  Poisson  distribution  is 
characterized  by: 

p(g,V  =Aa  'e  A/g[  [2] 
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For  large  A,  large  values  of  g  are  more  likely,  and  for  small  A,  small  values  of  g  are  more 
likely.  Therefore,  the  probability  of  g  can  be  modulated  by  varying  the  value  of  A 
according  to  a  sine  or  exponential  decay  function,  and  the  probability  of  small  gaps  can 
be  increased  where  the  SNR  is  highest  in  the  MR  signal  envelope  [32].  To  generate  g  as 
a  function  of  A,  a  Poisson  process  can  be  simulated  using  various  techniques  that  do  not 
depend  on  a  priori  knowledge  of  g  as  above  [37].  For  these  experiments,  the  poissrnd(A) 
function  in  Matlab  was  used  to  generate  g  as  a  function  of  A.  It  takes  as  input  an  array  of 
A  and  returns  an  array  of  sample  gaps  with  local  mean  and  variance,  A. 

An  example  of  the  2D  A,  gap,  and  mask  arrays  generated  by  Poisson-gap  mask  creation 
are  shown  in  Figure  2  and  illustrates  their  relationships;  as  the  size  of  A  and  the  gaps 
increase,  the  sample  density  decreases  in  that  area  of  the  mask.  The  magnitude  Point 
Spread  Function  (PSF)  of  each  mask  is  shown  and  demonstrates  the  viability  of  this 
approach;  the  single  dominant  central  peak  with  small  side-lobes,  surrounded  by  low 
amplitude,  incoherent  artifacts  is  the  desired  profile  of  a  random  under-sampling  mask 
PSF  [38]  [39], 

4D  MaxEnt  Reconstruction 

MaxEnt  image  reconstruction  was  used  in  these  experiments  to  reconstruct  the  under¬ 
sampled  data  sets.  MaxEnt  reconstruction  maximizes  the  entropy  of  the  data  in  the 
spatial,  spectral-domain  while  maintaining  data  fidelity  with  the  sampled  data  in  the  k- 
space,  time  domain  [25]. 

MaxEnt  image  reconstruction  of  4D  MRSI  data  is  formulated  as  a  constrained  convex 
optimization  problem  [25]  [27]: 

maximize  S(rri) iy  (3) 

s.  t  \\KTm  —  rf 111  <  C0 
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where  m  —  (y,  x,  F2,  F:)  is  the  reconstructed  spatial,  spectral-domain  data,  T  is  the  4D 
Fourier  operator,  K  is  the  under-sampling  mask  that  determines  which  samples  were 
acquired  in  the  ky-ti  plane,  d  =  (ky,  kx,  t2,  tx)  is  the  k-space,  time-domain  sampled  data, 
CQ  is  the  standard  deviation  of  the  noise  in  d,  and  S(m)  ij  is  the  entropy  of  the  estimated 
spectrum 


The  entropy  used  in  these  experiments  was  not  the  often  used  —  S  P  log(p)  entropy 
introduced  by  Shannon  [28]  but  the  Sx/2  entropy  derived  by  Daniell  and  Hore  specifically 
for  NMR  spectra  originating  from  spin  x/2  nuclei,  such  as  'fl  used  in  MRSI  [26]: 

( \m\ 


Z\m\ 

difloa 


7 


def 


+  4  + 


\m\j 


\\ 


def 


V 


—  4  + 


def , 


(4) 


) 


where  def  is  a  scaling  parameter  related  to  the  sensitivity  of  the  scanner  and  is  calculated 


for  m  of  length  N  as  j0/ ^  [27].  The  underlying  physical  processes  that  produce  an  MR 


spectrum  are  not  based  on  particle,  or  discrete  events  so  they  cannot  be  modeled  by 
simple  Poisson  processes  as  required  for  the  derivation  of  Shannon  entropy  [40],  They 
are  governed  by  the  density  matrix  of  the  spin  system  under  investigation.  This  equation 
can  be  applied  to  any  MR  spectrum  originating  from  spin  1/2  nuclei  and  addresses 
previous  concerns  regarding  the  use  of  entropy  in  MRS  and  MRI  reconstruction  [41], 


The  MaxEnt  reconstruction  problem  was  solved  by  a  Matlab  implementation  of  the 
Cambridge  algorithm  [25],  It  recasts  MaxEnt  image  reconstruction  into  an  unconstrained 
convex  optimization  problem  and  uses  a  variant  of  the  conjugate  gradient  method  to 
iteratively  find  the  extrema  in  two  phases;  the  first  phase  minimizes  the  fidelity 
constraint,  and  the  second  phase  maximizes  the  entropy,  while  keeping  the  fidelity 
constraint  minimized.  The  stopping  criterion  for  the  problem  is  reached  when  the 

gradients  of  the  entropy,  S(m),  and  fidelity  constraint,  C(m ),  are  parallel:  1^7 —  titH  < 


vs 

VC 

IISII2 

I|C||2 
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.001.  Specific  details  on  the  algorithm  and  modifications  to  accommodate 
multidimensional  MR  data  can  be  found  in  [27]. 

RESULTS 
MR  Simulations 

Quantitative  results  for  the  MaxEnt  reconstruction  of  the  simulated  4D  EP-COSI  data  set 
at  different  percent  under-samplings  and  SNR  are  shown  in  Figure  3.  As  described 
earlier,  noise  was  added  to  the  noiseless  4D  data  set  in  order  to  simulate  SNRs  of  2 
through  20  in  increments  of  2  which  were  then  undersampled  to  80%,  60%,  40%,  and 
20%  and  reconstructed  by  MaxEnt.  This  was  repeated  20  times  per  SNR  and  percent 
under-sampling  in  order  to  produce  a  mean  and  standard  deviation  for  each  metric. 

The  left  pane  of  Figure  3  shows  the  average  RMSE  versus  SNR  for  the  under-sampled 
and  MaxEnt  reconstructed  data  sets  at  each  percent  under-sampling.  The  RMSE  provides 
an  estimate  of  reconstruction  accuracy  with  respect  to  a  reference  data  set  that  increases 
as  the  two  data  sets  become  more  dissimilar.  The  RMSE  was  calculated  in  the  spatial- 
spectral  domain  as: 

RMSE  =  ^J^Qdataj  - \full\)2  (5) 

where  N  is  the  number  of  data  points,  full  is  the  fully  sampled  data  set,  and  data  is  the 
under-sampled  or  MaxEnt  reconstructed  data  set.  Error  bars  are  not  shown  because  the 
standard  deviations  were  three  to  four  orders  of  magnitude  smaller  than  the  mean  RMSEs 
and  did  not  vary  noticeably  over  percent  under-sampling  or  SNR.  As  can  be  seen,  the 
RMSE  of  the  under-sampled  data  set  increases  as  the  percent  of  sampled  data  decreases 
but  does  not  vary  considerably  with  SNR.  The  MaxEnt  reconstructions  decrease  the 
RMSE  significantly  at  each  SNR  and  percent  under-sampling  compared  to  the  under¬ 
sampled  data  set  but  begin  to  rise  at  very  low  SNR. 

The  right  pane  of  Figure  3  shows  the  average  NAA  to  Cho  and  Lac  to  Cr  ratios  versus 
SNR  for  the  fully  sampled,  under-sampled,  and  MaxEnt  reconstructed  data  sets  at  80%, 
60%,  40%,  and  20%  under-sampling.  The  ratios  were  calculated  from  the  peak  integrated 
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areas  of  each  metabolite  at  the  ppm  locations  listed  in  Table  1  and  then  normalized  by  the 
fully  sampled  ratios  at  SNR=20.  Each  peak  area  was  calculated  only  over  the  four 
spatially  distributed  voxels  for  each  metabolite.  Therefore,  the  NAA  peak  area  was 
calculated  over  voxels  (2,2),  (2,3),  (3,2)  and  (3,3)  as  shown  in  Figure  1,  for  example.  This 
ensured  that  any  spectral-spatial  leakage  caused  by  the  under-sampling  into  adjacent 
voxels  would  modify  the  peak  area  for  each  metabolite  and  change  the  metabolite  ratios 
as  their  peaks  aliased  along  the  ky-ti  plane. 

The  NAA/Cho  and  Lac/Cr  ratios  for  the  under-sampled  data  set  shown  in  the  bottom  of 
the  right  pane  in  Figure  3  differ  noticeably  from  the  fully  sampled  ratios  at  each  SNR. 

The  NAA/Cho  ratios  are  1-8%  larger  than  the  fully  sampled  ratios  and  increase  as  the 
percent  under-sampling  decreases.  The  Lac/Cr  ratios  decrease  as  less  data  is  sampled  and 
are  up  to  9%  lower  than  the  fully  sampled  ratios  for  sample  rates  of  80,  60,  and  40%,  and 
17%  lower  for  the  20%  sample  rate.  Both  the  NAA/Cho  and  Lac/Cr  ratios  show 
significant  deviations  in  the  metabolite  ratios  compared  to  the  fully  sampled  true  ratios; 
however,  these  deviations  are  limited  to  the  percent  under-sampling  and  are  not  largely 
affected  by  SNRs  down  to  2. 

The  MaxEnt  reconstruction  ratios  for  NAA/Cho  and  Lac/Cr  shown  in  the  top  of  the  right 
pane  in  Figure  3  show  much  better  agreement  with  the  fully  sampled  ratios  at  each  SNR 
and  percent  under-sampling  than  the  under-sampled  data.  The  NAA/Cho  ratios  are  within 
1%  of  the  fully  sampled  ratios  at  SNRs  greater  than  6  but  do  deviate  as  high  as  3%  for 
lower  SNRs.  The  Lac/Cr  ratios  are  also  improved  across  all  SNRs  to  be  within  5%  of  the 
fully  sampled  ratios  for  sample  rates  of  80,  60,  and  40%  and  12%  for  the  20%  sample 
rate;  however,  unlike  the  NAA/Cho  ratios,  the  Lac/Cr  ratios  were  affected  by  decreasing 
SNR  and  improvements  over  the  under-sampled  data  set  are  less  significant  at  SNRs  less 
than  6. 

Retrospective  Under-sampling  of  4D  EP-COSI 

The  results  from  retrospectively  under-sampling  to  25%  and  then  MaxEnt  reconstructing 
a  fully  sampled  4D  EP-COSI  scan  of  healthy  breast  in  vivo  are  shown  in  Figure  4.  The 
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mask  used  to  under-sample  the  ky-ti  plane  is  shown  on  the  far  left  along  with  the  signal 
envelopes  for  each  dimension  and  was  generated  using  the  2D  Poisson-gap  method 
described  earlier  in  this  paper.  The  central  point  in  ky  was  sampled  for  each  ti  increment 
and  the  1st  ti  point  was  sampled  at  each  ky.  [35]  [17].  Because  it  is  a  4D  data  set  where 
each  voxel  is  a  2D  spectrum,  space  limitations  prevent  us  from  showing  every  voxel  in 
the  16x16  FOV.  Therefore,  only  a  selected  2D  COSY  spectrum  from  the  data  set  as  well 
as  the  spatial  distribution  of  specific  spectral  peaks  is  shown. 

Figure  4  A1  shows  a  2D  COSY  spectrum  extracted  from  the  fully  sampled  4D  EP-COSI 
scan  of  healthy  breast  in  vivo.  It  was  taken  from  the  fatty  breast  region  highlighted  in 
Figure  4  A2.  It  clearly  shows  the  lipid  diagonal  peaks:  olefinic  fat  (UFD),  methyl  fat 
(FMETD),  and  fat  (FAT/FAT2/FAT3)  and  the  cross  peaks:  unsaturated  fatty  acid  right 
(UFR),  unsaturated  fatty  acid  left  (UFL),  and  triglyceryl  fat  (TGFR)  [7].  The  spatial 
distribution  of  the  UFL  and  UFR  cross-peaks  from  the  fully  sampled  4D  EP-COSI  scan  is 
shown  in  Figure  4  A2  overlaid  on  the  anatomical  MRI. 

Figure  4  B1  shows  the  same  2D  COSY  spectrum  as  Figure  4  A1  after  25%  under¬ 
sampling  has  been  applied  to  the  ky-fi  plane  using  the  mask  shown  at  left.  The  spatial- 
spectral  incoherent  artifacts  from  the  under-sampling  manifest  as  smeared  peaks  along  fi, 
which  is  illustrated  by  the  collapse  of  the  peaks  in  the  ID  projection  along  the  Fi 
dimension  on  the  right.  The  aliasing  of  the  large  diagonal  fat  peaks  around  (F2=2ppm, 
Fi=2ppm)  obscure  the  much  smaller  UFL  and  UFR  cross-peaks  around  (F2=2.1ppm, 
Fi=5.4ppm).  Figure  4  B2  shows  the  spatial  distribution  of  the  UFL  and  UFR  cross-peaks 
and  how  spatial  artifacts  from  under-sampling  ky-fi  manifest  as  errant  peaks  in  adjacent 
voxels  when  compared  to  the  distribution  of  fully  sampled  peaks  in  Figure  4  A2. 

Figure  4  Cl  shows  the  2D  COSY  spectrum  from  Figure  4  B1  after  MaxEnt 
reconstruction  has  been  applied.  The  majority  of  spectral  artifacts  seen  in  Figure  4  B1 
have  been  removed  and  the  metabolite  peaks  have  been  rebuilt  in  the  ID  projection  on 
the  right  of  the  figure.  The  aliasing  caused  by  the  diagonal  fat  peaks  around  (F2=2ppm, 
Fi=2ppm)  has  been  removed  and  the  UFL  and  UFR  cross-peaks  around  (F2=2.1ppm, 
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Fi=5.4ppm)  are  fully  resolved.  When  compared  to  the  fully  sampled  spectrum  in  Figure  4 
Al,  the  MaxEnt  reconstructed  spectrum  has  equivalent  spectral  peaks,  and  all  of  the 
significant  diagonals  and  cross-peaks  are  fully  resolved.  The  errant  peaks  in  the  spatial 
distribution  in  Figure  4  B2  caused  by  the  spatial  under-sampling  have  been  removed  in 
Figure  4  C2. 

The  range  of  values  for  the  (UFL+UFR)/(FAT3+FAT2)  integrated  peak  area  ratio  and 
average  error  from  the  fully  sampled  ratios  for  all  24  voxels  in  the  localized  VOI  is 
shown  in  Table  2  for  the  under-sampled  and  MaxEnt  reconstructed  data  sets.  As  can  be 
seen,  under-sampling  caused  the  ratios  to  vary  considerably  from  the  fully  sampled  ratios 
with  a  high  mean  ratio  error.  The  MaxEnt  reconstructed  ratios  show  much  better 
agreement  with  the  fully  sampled  ratios  and  have  a  much  smaller  mean  ratio  error. 

Prospectively  Undersampled  4D  EP-JRESI 

The  results  from  the  prospectively  under-sampled  4D  EP-JRESI  scan  of  a  human  prostate 
with  malignant  lesions  in  the  left  and  right  base  are  shown  in  Figures  5  and  6.  The  25% 
under-sampling  mask  with  signal  envelopes  along  the  ky  and  ti  dimensions  is  shown  on 
the  left  of  Figure  5.  Select  2D  JPRESS  spectra  from  the  4D  EP-JRESI  data  set  are  shown 
in  Figure  6  with  the  spatial  distribution  of  specific  spectral  peaks  in  Figure  5. 

Figure  5  shows  the  spatial  distributions  of  the  Cit  peaks  highlighted  in  Figure  6  for  the 
under-sampled  and  MaxEnt  reconstructed  data  sets;  the  purple  highlighted  voxel 
corresponds  to  the  location  of  the  extracted  2D  J-resolved  spectrum  in  Figure  6A,  the  red 
corresponds  to  Figure  6B,  and  the  blue  corresponds  to  Figure  6C.  Because  of  the 
sensitivity  profile  of  the  endorectal  coil,  the  SNR  of  voxels  close  to  the  rectal  wall  is 
higher  than  the  SNR  of  voxels  further  away  [42],  The  MaxEnt  reconstruction  in  Figure 
5B  shows  the  SNR  of  the  Cit  peaks  in  healthy  tissue  decreasing  further  from  the  rectal 
wall,  as  expected.  However,  the  spatial  distribution  of  the  under-sampled  data  set  in 
Figure  5A  shows  aliased  Cit  peaks  near  the  apex  of  the  prostate  as  well  as  high  SNR 
peaks  within  the  rectum  that  should  not  be  there.  The  spatial  distribution  of  Cit  in  the 


Magnetic  Resonance  in  Medicine 


Magnetic  Resonance  in  Medicine 


Page  14  of  30 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 


14 


under-sampled  data  is  noisier  in  general  and  shows  significant  spectral-spatial  artifacts 
when  compared  to  the  MaxEnt  reconstruction. 

The  top  of  Figure  6  shows  results  for  the  25%  under-sampled  4D  EP-JRESI  scan  and  the 
bottom  shows  results  for  the  MaxEnt  reconstruction  of  that  scan.  Because  the  scan  was 
prospectively  undersampled,  there  is  no  fully  sampled  data  for  comparison. 

Figure  6B  shows  2D  J-resolved  spectra  extracted  from  the  4D  NUS  EP-JRESI  data  sets 
that  are  highlighted  in  red  in  Figure  5.  As  shown  in  the  under-sampled  spectrum  in  the 
top  of  the  figure,  the  citrate  (Cit)  diagonal  and  its  cross  peaks  around  F2=2.6ppm  are  not 
easily  resolved.  There  is  significant  aliasing  of  the  Cit  peaks  along  fi  as  well  as  the  fat 
peaks  around  F2=1.7ppm.  The  bottom  of  the  figure  shows  the  same  2D  J-resolved 
spectrum  after  MaxEnt  reconstruction  of  the  4D  NUS  EP-JRESI  data  set.  The  Cit 
diagonal  and  cross  peaks  are  now  fully  resolved  with  no  aliasing  along  fi  as  well  as  the 
fat  peaks  around  F2=l  ,7ppm. 

A  high  Cit  to  Cho  ratio  is  nonnal  in  healthy  prostate  tissue,  but  a  low  Cit  to  Cho  ratio  is 
indicative  of  pathology  [43],  Therefore,  Figures  6A  and  6C  show  2D  J-resolved  spectra 
from  the  4D  EP-JRESI  data  set  from  tumors  in  the  left  and  right  base  of  the  prostate  that 
are  highlighted  in  purple  and  blue  in  Figure  5,  so  they  should  have  decreased  Cit  and 
elevated  Cho.  The  presence  of  tumors  was  verified  by  biopsy  with  Gleason  scores  of  3+3 
and  prostatic  specific  antigen  of  8.  The  under-sampled  spectra  in  the  top  of  Figure  6C 
show  Cho  peaks  at  F2=3.2ppm  that  are  highly  aliased  along  fi  and  Figure  6A  shows 
aliased  Cit  at  F2=2.6ppm  from  adjacent  voxels.  The  MaxEnt  reconstructions  are  shown  in 
the  bottom  of  the  figure  for  the  same  2D  J-resolved  spectra.  The  Cho  peaks  are  no  longer 
aliased  along  fi  and  the  errant  peaks  have  been  removed. 

DISCUSSION 

This  work  addresses  the  long  scan  time  requirements  of  4D  MRSI,  which  is  one  of  the 
main  issues  preventing  its  widespread  adoption.  Non-uniform  under-sampling  was  used 
to  accelerate  the  4D  data  acquisition  and  MaxEnt  non-linear  reconstruction  was  used  to 
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rebuild  the  missing  samples.  This  technique  resulted  in  minimal  losses  in  spatial-spectral 
resolution  while  reducing  scan  times  considerably. 

The  simulated  4D  EP-COSI  data  set  results  shown  in  Figure  3  demonstrate  how  the 
MaxEnt  reconstruction  perfonns  at  different  levels  of  noise  and  percent  under-sampling. 
The  MaxEnt  reconstruction  decreases  the  RMSE  significantly  at  each  SNR  and  percent 
under-sampling  compared  to  the  under-sampled  data  set  but  begins  to  rise  at  very  low 
SNR,  indicating  the  inability  of  MaxEnt  to  fully  reconstruct  the  data  set  when  features 
and  sampling  artifacts  are  obscured  by  high  levels  of  noise.  The  data  fidelity  constraint  in 
equation  (3)  determines  how  closely  the  MaxEnt  reconstructed  points  must  be  to  the 
sampled  points  to  within  the  standard  deviation  of  noise,  which  increases  as  the  noise 
floor  increases.  This  increase  in  the  noise  floor  effectively  “loosens”  the  fidelity 
constraint,  which  allows  the  reconstructed  points  to  vary  from  their  sampled  counterparts 
and  increases  the  entropy  of  the  reconstructed  spectrum.  As  the  entropy  increases,  so  too 
do  the  non-linearity  and  RMSE  of  the  reconstruction  with  a  loose  fidelity  constraint  [44], 
The  data  fidelity  constraint  can  be  “tightened”  beyond  the  standard  deviation  of  the  noise 
in  an  effort  to  reduce  the  RMSE  and  reconstruction  non-linearity  but  this  prevents 
MaxEnt  from  completely  removing  the  spatial-spectral  under-sampling  artifacts  close  to 
the  noise  floor,  which  could  potentially  obscure  small  features  [27]. 

The  results  from  Figure  3  show  the  NAA/Cho  ratios  are  consistent  with  the  fully  sampled 
data  set  ratios  and  offer  significant  improvement  over  the  under-sampled  data,  even 
without  using  methods  to  reduce  the  reconstruction  non-linearity  [45];  however,  the 
Lac/Cr  ratios  deviate  slightly  from  the  fully  sampled  ratios  at  higher  percent  under¬ 
sampling  and  significantly  at  20%  under-sampling.  Upon  further  investigation,  it  was 
found  that  each  metabolite  area  used  to  calculate  the  NAA/Cho  and  Lac/Cr  ratios,  with 
the  exception  of  Lac,  was  stable  and  consistent  with  the  fully  sampled  areas  over  SNR 
and  percent  sampling;  however,  the  areas  calculated  for  Lac  were  consistently  low  across 
all  SNRs  for  each  percent  under-sampling.  This  was  caused  by  a  mismatch  between  the 
2D  Poisson-gap  sample  masks  that  were  sine-bell  modulated  along  fi  and  the  filtered 
signal  envelope  of  Lac,  which  was  thought  to  be  a  sine-bell  but  upon  further  investigation 
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it  was  found  to  be  a  shifted  sine-bell.  This  shift  in  the  Lac  signal  envelope  meant  the 
mask  failed  to  sufficiently  sample  the  high  SNR  Lac  points  along  ti  for  each  percent 
under-sampling  and  prevented  the  full  Lac  peak  area  from  being  reconstructed  by 
MaxEnt  [45]  [35]  [34],  This  discrepancy  emphasizes  the  importance  of  the  sample  mask 
to  the  SNR  of  each  reconstructed  metabolite  for  4D  EP-COSI  data  and  the  dependence  on 
the  shape  of  standard  filters  applied  prior  to  reconstruction  [46].  For  each  metabolite  in 
EP-JRESI  data,  however,  the  maximum  SNR  point  along  fi  is  always  at  fi=l  and  is  not 
apodized  by  spectral  filters  prior  to  reconstruction;  therefore,  an  exponentially  damped 
sample  mask  that  only  considers  T2  signal  losses  in  the  time  domain  can  be  used  to 
sample  the  high  SNR  points  for  each  metabolite  without  the  need  to  consider  coherence 
transfer  signal  envelopes. 

The  Poisson-gap  sampling  masks  used  in  these  experiments  were  generated  by  a  random 
Poisson  distribution,  which  injects  a  degree  of  uncertainty  into  the  MaxEnt 
reconstruction.  It  has  proven  to  be  a  robust  technique  that  generates  masks  with  desirable 
PSFs,  as  shown  in  Figure  2,  and  the  RMSEs  of  reconstructed  data  sets  using  different 
Poisson-gap  sample  masks  are  stable  [32].  Previous  attempts  by  our  group  to  use 
detenninistic  masks  that  were  not  randomly  generated  were  difficult  to  optimize  and 
suffered  from  coherent  aliasing  which  cannot  be  removed  by  MaxEnt  reconstruction  [47] 
[35].  Recent  work  in  under-sampled  multidimensional  NMR  data  sets  using  detenninistic 
sample  masks  has  shown  promise  and  could  be  adapted  to  4D  MRSI  [48].  However, 
these  sample  masks  are  modulated  by  an  exponential  decay  function,  which  is  suitable  for 
EP-JRESI  data,  but  would  need  to  be  modified  to  use  a  sine-bell  function  for  EP-COSI 
data.  Poisson-gap  sampling  has  the  advantage  that  it  is  easily  modified  to  accommodate 
different  modulation  functions  and  can  therefore  be  used  for  many  different  scan  types.  It 
is  unknown  whether  these  new  detenninistic  methods  are  as  accommodating  as  Poisson- 
gap  in  that  regard. 

The  healthy  human  breast  results  in  Figure  4  show  that  under-sampling  and  MaxEnt 
reconstruction  work  well  for  filtered  in  vivo  EP-COSI  scans  with  a  sine-bell  signal 
envelope.  The  SNR-loss  caused  by  To  decay  that  was  not  present  in  the  simulated  4D  EP- 
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COSI  data  set  did  not  reduce  the  efficacy  of  the  reconstruction.  There  was  still  sufficient 
SNR  in  the  time  domain  to  reconstruct  the  diagonal  and  cross  peaks  in  the  spectral 
domain.  The  improvement  in  the  ratios  of  certain  lipid  peaks  shown  in  Table  2  indicates 
the  signal  from  the  large  diagonal  fat  peaks  that  had  been  aliased  over  the  ky-ti  plane  by 
the  under-sampling  was  returned  to  the  diagonal  peaks  by  MaxEnt  reconstruction  and 
surfaced  the  smaller  UFL  and  UFR  cross  peaks  that  had  been  obscured.  The  NUS  data  set 
in  Figures  4B1  and  4B2  show  artifacts  spread  along  fi  and  ky  as  well  as  reduced  spectral 
resolution  along  fi  caused  by  the  broad  NUS  PSF.  The  homogenous  nature  of  fatty  breast 
spectra  made  it  difficult  to  detennine  from  the  figure  if  the  spatial  resolution  along  ky 
decreased  as  a  result  of  the  NUS  PSF  as  expected.  However,  it  is  clear  the  effects  of  the 
NUS  PSF  along  fi  were  removed  by  MaxEnt  in  Figure  4C1  and  the  errant  spectral  peaks 
in  the  spatial  distribution  were  removed  in  Figure  4C2,  suggesting  the  spatial  PSF  along 
ky  was  improved.  Any  spectral  bleed  from  the  spatial  PSF  of  the  EP-COSI  pulse 
sequence  along  kx  was  orthogonal  to  the  effect  of  the  NUS  PSF  along  ky  and  was  not 
affected  by  the  MaxEnt  reconstruction. 

Under-sampling  of  4D  MRSI  data  can  be  implemented  in  the  pulse  sequence  on  the 
scanner  and  successfully  reconstructed  using  MaxEnt  as  shown  in  Figures  5  and  6  for  the 
prospectively  under-sampled  EP-JRESI  prostate  data.  Spatial-spectral  artifacts  caused  by 
the  broad  NUS  PSF  were  removed  by  MaxEnt  to  reveal  spectra  that  were  spatially 
localized  in  ky  properly  and  lacked  significant  spectral  aliasing  in  ti.  The  Cit  multiplet  in 
the  MaxEnt  reconstructed  data  set  showed  a  nicely  J-resolved  AB-type  structure  that  had 
been  aliased  by  convolution  with  the  NUS  PSF  in  the  under-sampled  data.  No 
comparisons  can  be  made  to  the  fully  sampled  data  set  as  was  done  for  the  retrospective 
EP-COSI  data  set;  however,  the  mean  (Cr+Cho+Spm)/Cit  ratio  for  the  MaxEnt 
reconstructed  2D  JPRESS  voxel  shown  in  Figure  6B  and  the  adjacent  four  voxels  within 
the  prostate,  excluding  the  tumors  in  the  left  and  right  base,  is  0.547  with  a  coefficient  of 
variation  of  21.2%  which  is  consistent  with  reported  literature  [43]  [49]. 

This  work  has  demonstrated  the  potential  for  MaxEnt  to  accelerate  the  acquisition  of  4D 
MRSI  data  in  vivo  into  clinically  acceptable  scan  times.  However,  because  of  its 
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similarity  to  CS  reconstruction,  quantitative  comparisons  between  these  techniques  need 
to  be  made  in  order  to  detennine  the  superior  method  for  reconstructing  4D  MRSI  data 
and  will  be  a  topic  of  future  research.  Additionally,  detennining  whether  Poisson-gap  or 
detenninistic  sample  masks  should  be  used  and  optimizing  the  modulation  functions  for 
4D  EP-JRESI  and  EP-COSI  data  sets  and  different  spectral  filters  will  be  addressed  in 
future  papers. 

CONCLUSIONS 

We  have  shown  that  MaxEnt  can  successfully  reconstruct  under-sampled  4D  MRSI  data 
by  reconstructing  a  mixed  domain  spectral-spatial  plane.  Simulated  4D  MRSI  data 
provided  a  quantitative  characterization  of  the  MaxEnt  reconstruction  at  different  percent 
under-sampling  and  SNR  and  it  was  shown  that  in  vivo  4D  EP-COSI  and  4D  EP-JRESI 
scans  could  be  reconstructed.  This  acceleration  translates  into  a  clinically  viable  5  minute 
EP-COSI  breast  scan  and  6  minute  EP-JRESI  prostate  scan. 
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Table  1:  PPM  locations  for  metabolite  peaks  used  in  ratio  calculations 


NAA 

Cho 

Cr 

Lac 

Diagonals  (ppm) 

(2.0, 2.0), 

(2.5, 2.5),  (4.3, 
4.3) 

(3.2, 3.2), 

(3.5, 3.5), 

(4.0, 4.0) 

(3-1, 3-1), 

(3.9, 3.9) 

(1-3, 1.3), 

(4.1, 4.1) 

Cross  Peaks  (ppm) 

(2.5, 4.3), 

(4.3, 2.5) 

(3.5, 4.0), 

(4.0, 3.5) 

N/A 

(1.3, 4.1), 

(4.1, 1.3) 

Table  2:  Range  and  mean  error  of  (UFL+UFR)/(FAT3+FAT2)  integrated  peak  area 
ratios  for  fully  sampled,  NUS,  and  MaxEnt  reconstruction. 


Ratio  Range 

Mean  Ratio  Error 

Full 

MaxEnt 

NUS 

0.1840-0.0634 
0.1604  -0.0550 
0.4262  -0.2180 

N/A 

0.0188  ±  .0048 
0.1737  ±  .0459 

Magnetic  Resonance  in  Medicine 


Magnetic  Resonance  in  Medicine 


Page  24  of  30 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 


24 


FIGURE  LEGENDS 

Figure  1:  Simulated  quad  phantom  illustration.  (Top)  Spatial  distribution  of  an  NAA,  Cr, 
Cho,  and  Lac  diagonal  peak  when  fully  sampled.  (Bottom)  Spatial  distribution  of  the 
same  NAA,  Cr,  Cho,  and  Lac  diagonal  peak  as  above  with  25%  under-sampling  in  the  ky- 
ti  plane. 

Figure  2:  Poisson  sample  mask  creation  along  the  ky-ti  plane  for  EP-JRESI  and  EP- 
COSI.  (Top)  Sine  and  exponentially  modulated  values  for  A.  (Upper  Middle)  Poisson 
distributed  values  for  each  A  indicating  the  sample  gap  between  that  point  and  other 
points  in  the  sample  distribution.  (Lower  Middle)  Resulting  2D  Poisson-gapped  sample 
mask.  (Bottom)  Magnitude  PSF  of  the  2D  sample  mask. 

Figure  3:  Metrics  comparing  the  under-sampled,  MaxEnt  reconstruction,  and  fully 
sampled  4D  EP-COSI  simulated  data.  (Left)  RMSE  of  under-sampled  and  MaxEnt 
reconstructed  data  versus  SNR  for  20,  40,  60,  and  80%  under-sampling.  (Right) 
Comparison  of  the  NAA/Cho  and  Cr/Lac  ratios  for  fully  sampled  data  to  MaxEnt  (Top) 
and  under-sampled  (bottom)  reconstruction  versus  SNR  for  20,  40,  60,  and  80%  under¬ 
sampling 

Figure  4:  (Left)  Under-sampling  mask  used  to  retrospectively  under-sample  the  ky-ti 
plane  to  25%  (Top)  Selected  2D  COSY  spectrum  from  a  4D  EP-COSI  scan  of  healthy, 
fatty  breast  highlighted  in  the  bottom  images  for  Al)  fully  sampled  Bl)  25%  under¬ 
sampling  Cl)  MaxEnt  reconstruction  (Bottom)  Spatial  distribution  of  the  olefinic  cross¬ 
peaks  of  unsaturated  fatty  acids  (UFL/UFR)  highlighted  in  the  2D  COSY  spectrum  for 
A2)  fully  sampled  B2)  25%  under-sampling  C2)  MaxEnt  reconstruction. 

Figure  5:  Spatial  distribution  of  the  citrate  peaks  from  a  4D  EP-JRESI  scan  of  malignant 
prostate  highlighted  in  the  2D  JPRESS  spectra  from  Figure  6  for  (A)  25%  Under-sampled 
(B)  MaxEnt  reconstructed.  (Left)  Mask  used  to  prospectively  under-sample  the  ky-ti 
plane  to  25%. 

Figure  6:  Select  2D  J-resolved  spectra  from  a  4D  EP-JRESI  scan  of  malignant  prostate 
for  (top)  prospectively  25%  under-sampled  and  (bottom)  MaxEnt  reconstructed.  A) 
Spectra  highlighted  in  purple  in  Figure  5  showing  elevated  Choline  and  depressed  Citrate 
in  a  malignant  lesion.  B)  Spectra  highlighted  in  red  in  Figure  5  showing  elevated  Citrate 
and  nonnal  Choline  in  healthy  prostate.  C)  Spectra  highlighted  in  blue  in  Figure  5 
showing  elevated  Choline  and  depressed  Citrate  in  a  malignant  lesion. 
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Figure  1:  Simulated  quad  phantom  illustration.  (Top)  Spatial  distribution  of  an  NAA,  Cr,  Cho,  and  Lac 
diagonal  peak  when  fully  sampled.  (Bottom)  Spatial  distribution  of  the  same  NAA,  Cr,  Cho,  and  Lac  diagonal 
peak  as  above  with  25%  under-sampling  in  the  ky-tl  plane. 
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Figure  2:  Poisson  sample  mask  creation  along  the  ky-tl  plane  for  EP-JRESI  and  EP-COSI.  (Top)  Sine  and 
exponentially  modulated  values  for  A.  (Upper  Middle)  Poisson  distributed  values  for  each  A  indicating  the 
sample  gap  between  that  point  and  other  points  in  the  sample  distribution.  (Lower  Middle)  Resulting  2D 
Poisson-gapped  sample  mask.  (Bottom)  Magnitude  PSF  of  the  2D  sample  mask. 
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Figure  3:  Metrics  comparing  the  under-sampled,  MaxEnt  reconstruction,  and  fully  sampled  4D  EP-COSI 
simulated  data.  (Left)  RMSE  of  under-sampled  and  MaxEnt  reconstructed  data  versus  SNR  for  20,  40,  60, 
and  80%  under-sampling.  (Right)  Comparison  of  the  NAA/Cho  and  Cr/Lac  ratios  for  fully  sampled  data  to 
MaxEnt  (Top)  and  under-sampled  (bottom)  reconstruction  versus  SNR  for  20,  40,  60,  and  80%  under¬ 
sampling 
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Figure  4:  (Left)  Under-sampling  mask  used  to  retrospectively  under-sample  the  ky-tl  plane  to  25%  (Top) 
Selected  2D  COSY  spectrum  from  a  4D  EP-COSI  scan  of  healthy,  fatty  breast  highlighted  in  the  bottom 
images  for  Al)  fully  sampled  Bl)  25%  under-sampling  Cl)  MaxEnt  reconstruction  (Bottom)  Spatial 
distribution  of  the  olefinic  cross-peaks  of  unsaturated  fatty  acids  (UFL/UFR)  highlighted  in  the  2D  COSY 
spectrum  for  A2)  fully  sampled  B2)  25%  under-sampling  C2)  MaxEnt  reconstruction. 
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Figure  5:  Spatial  distribution  of  the  citrate  peaks  from  a  4D  EP-JRESI  scan  of  malignant  prostate  highlighted 
in  the  2D  JPRESS  spectra  from  Figure  6  for  (A)  25%  Under-sampled  (B)  MaxEnt  reconstructed.  (Left)  Mask 
used  to  prospectively  under-sample  the  ky-tl  plane  to  25%. 
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Figure  6:  Select  2D  J-resolved  spectra  from  a  4D  EP-JRESI  scan  of  malignant  prostate  for  (top)  prospectively 
25%  under-sampled  and  (bottom)  MaxEnt  reconstructed.  A)  Spectra  highlighted  in  purple  in  Figure  5 
showing  elevated  Choline  and  depressed  Citrate  in  a  malignant  lesion.  B)  Spectra  highlighted  in  red  in  Figure 
5  showing  elevated  Citrate  and  normal  Choline  in  healthy  prostate.  C)  Spectra  highlighted  in  blue  in  Figure  5 
showing  elevated  Choline  and  depressed  Citrate  in  a  malignant  lesion. 
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